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Theoretical antiproton and proton cross sections for ionization and excitation of hydrogen 
molecules as well as energy spectra of the ionized electrons were calculated in the impact-energy 
range from 8 to 4000 keV. The cross sections were computed with the close-coupling formulation of 
the semi-classical impact-parameter method. The target was described using a one-active electron 
model centered on the target and assuming a fixed internuclear distance during the collision pro- 
cess. The dependence of the ionization cross sections on the internuclear distance is examined. The 
present cross sections are compared with experimental and theoretical data from the literature. For 
impact energies _E > 90 keV the obtained results for ionization by antiproton impact are comparable 
to the experimental data while they disagree for energies -E < 80 keV. 



PACS numbers: 34.50.-s,34.50.Gb 



I. INTRODUCTION 



The ionization of atoms and molecules by antiproton 
[p) impacts has become the subject of great theoreti- 
cal and experimental interest. The motivation for this 
is twofold. First, there is a fundamental interest in colli- 
sions involving exotic particles and second, in comparison 
to proton {p) collisions the sign of the projectile charge 
is exchanged which opens up the possibility to explore 
interesting physical effects. Further attention is drawn 
to this topic due to the upcoming Facility for Antiproton 
and Ion Research (FAIR) ^]. The international collabo- 
rations on atomic and molecular physics FLAIR [2] and 
SPARC [3], both being a part of the FAIR project, intend 
to investigate antiproton driven ionization processes and 
even kinematically complete antiproton collision exper- 
iments [^. These experimental efforts complement the 
recent intense activity in experimental [5J,|6| and theoret- 
ical 0, B i, 113, [m, [11 anti-hydrogen studies. 

In the last 15 years a large number of elaborate cal- 
culations have been carried out for ionization in the 
sim ples t one-electron and two-electron swtems p -f H 
H, M 111 [MIMIlI,[ll:,|2Ml|j2iJMJ3 and p -h He 



23,Ml2al30,l3lin3ll3ln3il3^ respectively, 
and have offered the single-ionization cross sections in 
agreement with experimental results [3a, [33, [sl, [3^, So, 
l4l| for incident energies E > 50keV. An interesting dif- 
ference between the proton and antiproton impacts was 
experimentally recognized for the double-ionization cross 
sections of the He atom [3g, |33, |4^ , and its origin could 
be explained by theoretical studies [3^, \^^ ■ At low en- 
ergies (E < 50keV), however, large discrepancies still 
remained between the experiment and theoretical results 
both for the single- and double-ionization processes. Very 
recent experimental results for p -\- He collisions at en- 
ergies below 30keV could partly resolve this discrepan- 
cies [431 . Theoretical work has also been done for other 
atomic targets, e.g., recently for alkali- metal atoms |44l |. 

For molecular targets measurements were further made 
for ionization in p + H2 collisions [40, |4^. Again, a no- 



table difference between proton and antiproton impacts 
could be seen for the ionization cross sections. While de- 
tailed work was done for proton impacts [4^, [43, [4^, [40, 
[50 | , little is investigated about the ionization of molecules 
by antiproton impact. In the case of p -I- H2 colli- 
sions a calculation was done by Ermolaev [5lj , who used 
an atomic hydrogen target with a scaled nuclear charge 
Zn = l-OOZproton in Order to mimic the H2 ionization 
potential. A molecular approach was used by Sakimoto 
for calculations of the one-electron system p -\- H J [5^ • 
The aim of the present work is to examine the p -\- 
H2 collision process in some detail and to improve the 
existing theoretical cross sections. Therefore, the target 
molecule is described with an improved model potential 
in comparison to the simple model used by Ermolaev 
in J5ll |. The present model provides cross sections for 
different internuclear distances of the H2 molecule. Fur- 
thermore, the description of the continuum is improved 
by expanding the wave functions in B-splinc functions as 
well as by increasing the number of basis states. The 
paper is organized as follows: Sec. [TT| considers the de- 
scription of the hydrogen molecule as well as the com- 
putational approach. Sec. IIIII reports on the dependence 
of the ionization cross sections on the molecular inter- 
nuclear distance. Subsequently, the calculated cross sec- 
tions for proton and antiproton impact are presented and 
compared to literature data. Finally, the energy spectra 
of the ejected electrons are taken into account. Sec. IIVI 
concludes on the present findings and discusses the appli- 
cability of the used model. Atomic units are used unless 
it is otherwise stated. 



II. METHOD 
A. Model potential 

The target molecule is treated as an effective one- 
electron system. The effective electron is exposed to a 



model potential 



Vmod{r) = 1 1 + 1— r exp 

r V \a\ 



2r 



11/2 



(1) 



which contains one dimensionless parameter a and the 
nuclear charge is Z„ = +1. It was proposed by Vanne 
and Saenz in [5J] and farther discussed in [53|- They 
successfully used Vmod{r) for calculating the ionization 
and excitation of H2 molecules in intense ultrashort laser 
pulses. A comparison with their results from a full molec- 
ular treatment of H2 confirmed the applicability of the 
model for ionization as well as for excitation. The model 
potential reduces to the ionization potential of atomic 
hydrogen, /p[H] = 0.5 a. u., for a — » 0. Furthermore, it 
satisfies the conditions Knod(?") -^ — l/r for r ^ 00 and 
therefore describes, in contrast to the model used by Er- 
molaev in |51| , the long-range behavior of the potential 
correctly. 

The dependence of the ionization potential /mod(c^) on 
the parameter a for a g^^tem described by Knod is deter- 
mined numerically in [SJ] but for < a < 0.3 it can also 
be approximated well with the analytic expression 



In 



(2) 



a is chosen such that /mod corresponds to the ionization 
potential /p[H2] of the H2 molecule at a given intcrnu- 
clear distance it!„. /p[H2] as a function of i?„ is obtained 
by subtracting the ground-state potential-energy curve of 
H2 which was very accurately calculated by Wolniewicz 
[53 | from the ground-state energies of H^. The values 
of a used in this work for the various -R„ ranging from 
1.0 a. u. to 2.11 a. u. are given in TableUtogether with the 
corresponding ionization potentials. It should be men- 
tioned that within this model the effect of anisotropy 
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Table I: Values of the model potential parameter a used in 
this work for different internuclear distances Rn given in a. u. 
For these internuclear distances also the ionization potential 
/p[H2] for H2 is given in Hartree which is obtained using the 
H2 ground-state potential-energy curve calculated by Wol- 
niewicz [5311. 



due to the two nuclei as well as of the second electron is 
solely contained as a screening of the Coulomb potential. 
A different model potential for the target can be ob- 
tained if in Eq. ([T]) the parameter a is set to a — > and 
the nuclear charge is scaled to Z„ = Zsca.\ — 1.09 as it 
was, for example, proposed by Ermolaev in [5l|. In what 
follows, it shall be referred to this special case of Vmod as 
14cai which describes a scaled hydrogen atom Hgcai with 
all energy levels shifted according to 



Cj [Hscal] = (-^scal) Ej [H] . 



(3) 



The ionization potential of Hgcai is equal to the absolute 
value of the ground-state energy Ip\Rsca\] = I ei Pscai] I = 
0.59405 a. u. and corresponds to the ionization potential 
of a H2 molecule at an internuclear separation Rn = 
1.474 a. u. 



B. Time propagation 

An approach similar to the one used in this work was 
already applied and described in some detail in a previ- 
ous work |44| dealing with slow antiproton and proton 
collisions with alkali-metal atoms. The relative motion 
of the heavy projectiles is approximated by classical tra- 
jectories with constant velocities Vp parallel to the z axis 
also known as the impact-parameter method. The dis- 
tance vector R between the projectile and the target is 
given by 



R(t) = b + vpi. 



(4) 



where b is the impact-parameter vector along the x axis 
and t the time. The internuclear distance i?„ of the two 
molecular target nuclei is held fixed during the collision 
process. The time-dependent Schrodinger equation 



Vi„t(r,R(t)) vl/(r,i) (5) 



^^*(r,i)=(i^o 



of the target interacting with the projectile is solved. Hq 
is the Hamiltonian of the effective target electron includ- 
ing the model potential (H]) 



Ho = 



1. 



Vr 



mod 



(6) 



The time-dependent interaction between the projectile 
with the charge Zp and the target is given by the time- 
dependent interaction potential 



Mnt(r,R(i)) = 



-Z„ 



Zn 



|r-R(t)| |R(t)| 



(7) 



where r is the spatial coordinate of the effective electron. 
The time-dependent wavefunction 



*(r,t) 



^ Cj(i)0j(r) exp[-ieji] (8) 



is expanded in eigenstates (j)j of the target Hamiltonian 
Hq with the energies Cj . The radial part of (j>j is further- 
more expanded in B-spline functions whereas its angular 
part is expanded in spherical harmonics. Substitution of 
the wavefunction ^ in Eq. ([5]) by its expansion given in 
Eq. ([U results for a trajectory with fixed b (of. Eq. (JU) 
in a system of coupled differential equations 

i—Cj{t,b) = ^Cfc(i,6)((/)j|V-nt|0fc)exp[i(ej-efe)i]. (9) 

k 

The differential equations ([9|) are solved with the initial 
conditions Cj{t — —oo;b) = Sji, where ci is the coeffi- 
cient of the ground state of the target Hamiltonian Hq. 
The transition probability Pj{b) for a transition into the 
target state j after a collision with the impact parameter 
b is given by 

P,{b) = \c,it = +^,b)f (10) 

and accordingly the probabilities for ionization Pionib) 
and excitation Pcx{b) by 



Pionib) 



PUb) 



ej>0 



\cj{t^ +oo,b)\' 



J2 \cjit = +cx,,b)\' 

ei<ej<0 



(11) 

(12) 



The cross section for the transition into the target state 
j can be expressed as 



CTj =2 7r f Pj{b)bdb. 



(13) 



For an H2 target the cross section for ionization is defined 
as 



fTio„ = 2^ / 2Pion(&)6dfe 



and for excitation as 



.„ = 2,/2P„(i).d.. 



The electron-energy spectrum 

dcTion 



S{e)^ 



£,-=£ 



piej)aj, 



(14) 



(15) 



(16) 



that is the differential cross section for ejecting an elec- 
tron with an energy e, can be obtained by using the 
density of the continuum states p. The factor two in 
Eqs. (fHHTB)) accounts for the two indistinguishable elec- 
trons of the H2 molecule. The introduction of such a 
simple factor is believed to work as long as indepen- 
dent interactions of the electrons with the projectile are 
dominating. This is the case for high projectile veloc- 
ities. However, it is unclear to what extend phenom- 
ena including electron-electron-interaction like double- 
ionization can be described with an effective one-electron 
model and furthermore can be extracted from the deter- 
mined transition probabilities (Eqs. (fTO|) to (fTT|) ) and 
therfore are not considered here. 
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Figure 1: Ionization cross section (Tion(-Rn) for p + H2 as a 
function of the internuclear distance i?„ at E — 16, 64, 125, 
250, 1000, and 4000 keV. 



III. RESULTS AND DISCUSSION 

A. Dependence of the internuclear distance 

The calculated cross sections depend primarily on the 
ionization potential which in turn depends for the H2 
molecule on the internuclear distance i?„ (cf. Table [l|. A 
cross section <7{Rn) which depends on i?„ may therefore 
be defined as 



a(i?„) = a(/p[H2](i?„)) 



(17) 



As has been shown for p + H^ collisions by Sakimoto 
[53 | who used a full description of the molecular tar- 
get ion the dependence of crion{Rn) on _R„ in the range 
1.5 a. u. < Rn < 3.0 a. u. differs for parallel and perpen- 
dicular orientations. For an orientation of the molecu- 
lar ion perpendicular to the collision plane the depen- 
dence on Rn is rather weak. Whereas for an orienta- 
tion parallel to the trajectory of the projectile the ion- 
ization cross sections increase approximately linearly for 
larger Rn by more than a factor two in the given i?„ 
range. A stronger orientational dependence is expected 
for HJ cations compared to a H2 molecule because of 
their larger equilibrium internuclear distance. Necessar- 
ily, no molecular-orientation dependence is taken into ac- 
count in the present study. Since the present cross sec- 
tions depend according to Eq.[T7]on the ionization poten- 
tial which is orientation-independent the present results 
are interpreted as orientation-averaged. 

Fig. [T] shows the present p -|- H2 ionization cross sec- 
tions aion{Rn) for the impact energies E = 16, 64, 125, 
250, 1000, and 4000 keV as a function of i?„. Ionization 
cross sections were calculated for eleven different internu- 
clear distances in the range 1.0 a. u. < i?„ < 2.11 a. u. It 



can be seen that for all impact energies shown here the 
ionization cross sections increase with larger i?„. This 
can be explained with the decrease of the ionization po- 
tential for an increasing internuclear distance also shown 
in Table IH The dependence of aion on _R„, however, 
diminishes with higher impact energies. For energies 
E > 1000 keV crion(^n) depends only weakly on i?„ and 
increases by about a factor 1.2 in the whole i?„ range. For 
smaller energies E < 125 keV the dependence on the in- 
ternuclear distance is much stronger and for E = 16keV 
crion(-Rn) increases by more than a factor 1.7 in the con- 
sidered Rn range. Furthermore, for all considered impact 
energies the cross sections show roughly a linear depen- 
dence on Rn. Therefore, one may assume that for all 
impact energies E considered in the present investigation 
the simple relation 
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fionC-Rn) = CTioniRn ) + {Rn — Rn 



daion{Rn) 



dRn 



(18) 



holds approximately in the examined interval where i?„ 
is a fixed internuclear distance within this interval. 

In Fig. [2] ionization probabilities Pionib) weighted with 
the impact parameter 6 as a function of b are com- 
pared for four different internuclear distances i?„ = 1.2, 
1.5, 1.8, and 2.11a. u. at two antiproton impact energies 
E—125 keV and 1000 keV. It can be seen that the curves 
for the higher impact energy _E=1000keV differ much 
less than those for i? = 125keV in accordance with Fig. 
[TJ All maxima of the curves for E ~ 1000 kcV lie around 
b — 1.3 a. u. The maxima for E = 125 keV slightly move 
from b — 1.0 a. u. for i?„ = 1.2 a. u. towards b — 1.3 a. u. 
for Rn = 2.11a. u. and thereby also increase in height. 
Whereas, the qualitative behavior of the curves for a con- 
sidered impact energy does not change for varying i?„. 

In order to determine results which include the rovi- 
brational motion of the H2 molecule one may use the fact 
that the cross sections to a given electronic state can be 
correctly given using closure (cf. [56[) 



a{Rn) 



S,o{Rn) 
Rn 



{Rn) dRn 



(19) 



where S,o{Rn)/Rn is the radial nuclear wave function of 
an II2 molecule in its rovibronic ground state. Clearly, 
the integration in Eq. p9|) leads to a loss of the electron- 
energy resolution. The energy information, however, is 
not relevant for integrated cross sections a but for dif- 
ferential cross sections like the electron-energy spectrum 
do-/de. 

It is always possible to express a{Rn) in terms of an 
(infinite) polynomial in i?„ and therefore to reformulate 



Figure 2: (Color online) Ionization probability Pion{b) 
weighted with the impact parameter b as a function of b 
at two different antiproton impact energies _E = 125keV and 
1000 keV. Curves for four different internuclear distances Rn 
are compared [in a.u.]: 1.2 (dot-dashed curve), 1.5 (solid 
curve), 1.8 (dashed curve), and 2.11 (dotted curve). 



Eq. ([T9D as 






fe=0 



fc=l 



Co(-Rn) 
flfc / (i?n)'' \URn)\^ 

ak { {Rn) '^ ) + ao , 



di?. 



dRn (20) 
(21) 
(22) 



where (^ {Rn) ^ ) denotes the expectation value of (i?„) '"' 
for the rovibrational ground state of H2. If the cross 
section a{Rn) depends on i?„ linearly, which is here at 
least to a good extend the case, one finds, using Eqs. (fT5)) 
and (|22p . the special relation 



ai (i?„) +ao = o-((i?„)) 



(23) 



i.e., it is sufficient to evaluate the cross section at the 
expectation value of the internuclear distance ( i?„ ) of 
the II2 molecule. The value (i?„) = 1.448 a. u. has been 
reported by Kolos and Wolniewicz [53| and it was used in 
the present calculations to determine the ionization and 
excitation cross sections. 

It may be mentioned that although the vibration and 
rotation of the H2 molecule is taken into account a distor- 
tion of the molecular vibration and rotation during the 
collision with the projectile may possibly lead to a sub- 
stantial change in magnitude of the cross section. The 
effect of such a distortion (which is not accounted for 
in the present work) on a may be largest for small im- 
pact energies where the cross section depends consider- 
ably on Rn as has been shown in Fig.[T] In order to better 
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Figure 3: (Color online) Electron-loss cross sections aioss for p 
+ H2 as a function of the impact energy E. Theory. Present 
results: solid curve, p + H2; dotted curve, p + H multi- 
plied by two; dash-dotted curve, p + Hscai multiplied by 
two. Elizaga et al. [50]: dashed curve, optimized dynamical 
pseudostates method; dash-double-dotted curve, molecular 
orbitals method; double-dash-dotted curve, eikonal classical 
trajectory Monte Carlo method. Experimental results: filled 
circles, Rudd et al. [48l.l49|: filled squares, Shah and Gilbody 

Mm- 



understand collision processes involving slow antiprotons 
{E < 100 keV) it would be desirable to fully include, in 
an advanced approach, the evolution of the internuclear 
distance during the collision. 



B. Ionization of H2 by p impact 

As has been discussed in a previous work [4j in some 
detail much more effort is needed to bring proton com- 
pared to antiproton cross sections to convergence us- 
ing the present method. This is in particular true for 
low proton impact energies where electron capture be- 
comes the dominant loss channel of the target electrons. 
The difficulties in the description of the electron capture 
are mainly due to the use of a one-center expansion of 
the basis around the target which has to be compen- 
sated with an enlarged basis set. The main motivation 
for the present calculations of proton results is given by 
the need for a comparison of the employed method with 
an extended amount of literature since the experimen- 
tal and theoretical data on antiproton collisions with H2 
molecules are still sparse. A one-center expansion around 
the target, however, seems to be justified for antiproton 
collisions in which electron capture is absent and which 
are in the focus of this investigation. 

The present results for the electron loss of molecular 



hydrogen in collisions with protons are shown in Fig. [3] 
as a solid curve. Also shown are the electron-loss cross 
sections for atomic hydrogen in a p + H collision multi- 
plied by two. The present data are compared with ex- 
perimental results by Rudd et al. [48,, ,42J and by Shah 
and Gilbody [Mlil- 

The present findings for H2 match the experimental 
data by Rudd et al. in the whole energy range. The 
agreement with the measurements of Shah and Gilbody 
is also good except for £■ < 20 keV where their data starts 
to be smaller than the results of the present work as well 
as those of Rudd et al. The electron loss cross sections for 
an atomic hydrogen target in p -|- H collisions multiplied 
by two agree well with the experimental and present data 
for E > 300 keV. With decreasing impact energies, i.e., 
with increasing dependence of the cross sections on the 
internuclear distance, the results for p + H get, however, 
considerably too large. 

In the theoretical work by Elizaga et al. [50| a simi- 
lar model potential was used which can be obtained by 
integrating an effective hydrogen atom-like charge distri- 
bution with Gauss's theorem. This model potential was 
also proposed by Hartree in [58| for He atoms (i?„ = 0). 
Cross sections for the electron loss were calculated for 
Rn ~ 1.4 a. u. Thereby, the three methods molecular or- 
bitals (MO), optimized dynamical pseudostates (ODP), 
and eikonal classical trajectory Monte Carlo (ECTMC) 
were used in the calculations and the results are also 
shown in Fig. [31 The cross sections obtained with ODP 
are very similar to the present ones. Only for E < lOkeV 
they are larger than the present data and those by Rudd 
et al. The MO approach was applied only at low ener- 
gies E < 25 keV and leads throughout to similar, though, 
slightly larger results than those obtained with ODP. Ex- 
actly in the latter energy range the outcome obtained 
with ECTMC differs considerably from all other curves 
whereas for E > 25 keV it matches the experimental and 
the present results very well. It can be concluded that 
the present approach is capable of describing collisions 
with II2 targets quite accurately in the considered en- 
ergy range. 



C. Ionization of H2 by p impact 

The present results for the ionization of molecular hy- 
drogen by antiproton impact are shown in Fig. [4| as solid 
curve. Also shown are the ionization cross sections for 
antiproton collisions with atomic hydrogen multiplied by 
two. The results are compared with calculations by Er- 
molaev [5l| and experimental data for non-dissociative 
ionization by Anderson et al. :38] as well as data of a sub- 
sequent measurement by Hvelplund et al. |40|. As has 
been suggested by Hvelplund et al. in [40| the data for 
impact energies below 200 keV of their earlier measure- 
ment nil are omitted in Fig.[ll The data for E < 200 keV 
of their first experiment are generally some 10% larger 
than those in f4(| but have a considerably lower accu- 
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Figure 4: (Color online) Ionization cross section aion for p 
+ H2 as a function of the impact energy E. Theory: solid 
curve, present results; dotted curve, present results for p + H 
multiplied by two; dash-dotted curve, present results for p + 
Hscai multiplied by two; dashed curve, results for p + Hscai 
multiplied by two by Ermolaev [5l|]. Experiment: circles, 
Anderson et al. [i^]; squares, Hvelplund et al. |40| ]. 
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Figure 5: (Color online) Excitation cross section CTcx for p + 
H2 as a function of the impact energy E. Theory: solid curve, 
present results; long-dashed curve, results for p + H; dotted 
curve, same results for p + H multiplied by two; dash-dotted 
curve, results for p + Hscai multiplied by two. Cross sections 
for excitation into I — 1 states. Theory: thin short-dashed 
curve, present results; thin dash-double-dotted curve, results 
for p + Hscai multiplied by two. 



racy. 

For liigh impact energies E > 1000 IceV all theoreti- 
cal curves coincide and also agree with the experimental 
data. For lower energies (400 keV < E < 1000 keV) the 
ionization cross sections for atomic hydrogen start to dif- 
fer from both theoretical results for a hydrogen molecule. 
However, at these energies the atomic results seem to de- 
scribe better the experimental data. In the energy regime 
from 250 keV down to 90keV the theoretical cross sec- 
tions by Ermolaev approach those of the p + il calcu- 
lation which differ significantly from the measured cross 
sections. The experimental data are, however, well de- 
scribed by the present p + H2 cross section in this energy 
regime. Though, the strong variation of the experimen- 
tal data around 85keV is not followed by the smooth 
curve of the present results. While the magnitude of 
the present cross sections is comparable to the experi- 
mental data down to 20keV the functional behavior of 
both, experimental and present curve, starts to differ for 
E < 50keV. Here, the present p -I- H2 curve possesses 
a similar characteristic as two times the cross sections 
of the hydrogen atom but with smaller magnitude be- 
cause of the larger ionization potential of the molecule. 
The slope of the present cross sections at these low ener- 
gies may indicate the lack of two-electron effects in the 
target description. The experimental data, on the other 
hand, show a behavior very similar to that of the single 
ionization of helium also measured with the same experi- 
mental set-up by Hvelplund et al. [4Q| . Very recently the 
same authors published another measurement of the sin- 
gle ionization cross section for p -t- H2 in the energy range 
3 keV < E < 25 keV 43,1 which revealed that their helium 



single ionization cross sections in [40| are too small for 
the lowest measured energies. It may be an interesting 
question whether the same is true in the case of the p + 
H2 ionization cross sections as suggested by the present 
results. Therefore, it would be worthwhile to initiate a 
further attempt to measure p -I- H2 cross sections at low 
antiproton energies. 

An effective one-electron description with a fixed in- 
ternuclear distance seem to be sufficient to describe non- 
dissociative ionization cross sections for p -j- H2 at high 
energies. But it is unclear how strong the infiuence of 
two-electron effects and the variation of the internuclear 
distance is at intermediate and low energies. Since the 
energy regime around and below the maximum of the 
ionization cross section is considered to contain interest- 
ing physical effects a full quantum mechanical treatment 
of the target molecule would be desirable. It should be 
mentioned, however, that such an approach is very de- 
manding. 



D. Excitation of H2 by p impact 

The present excitation cross sections for p -f H2 are 
shown in Fig. [5] as solid curve. Also shown are results for 
antiproton collisions with atomic hydrogen and the same 
atomic cross sections multiplied by two. To the best of 
the authors' knowledge there are no data in literature to 
compare these results with. 

Due to the experiences with the ionization cross sec- 
tions one may estimate the range of validity of the excita- 



tion cross sections presented here to be about 100 keV < 
E < 4000 keV. Comparing the results for ionization and 
excitation in p + H2 cohisions one can say that cr^x is 
smaller than CTion for impact energies E < 1000 keV and 
that both are practically the same for larger energies. 
The maximum of (7cx.{E) lies around E = 58keV and 
therefore at a higher energy than the maximum for ion- 
ization. 

The excitation cross sections for molecular hydrogen 
can also be compared with the results for atomic hydro- 
gen. Fig. [S] clearly shows that the naive assumption an 
H2 molecule is essentially composed of two independent 
hydrogen atoms yields excitation cross sections which are 
obviously different from those which were obtained with 
the model potential Vmod given in Eq. ([1]). Only for high 
impact energies both curves get close to each other. On 
the other hand, it is interesting to observe that the ex- 
citation cross sections for a single hydrogen atom seem 
to be much more in accordance with the present molec- 
ular (Tex- Both cross sections show the same behavior 
and have practically the same values in the considered 
energy range. This similarity for atomic and molecular 
hydrogen targets was evidently not found in the case of 
ionization in Sec. IIII CI 



E. Electron spectra 

The electron-energy spectra S{e, E) = dcr(e, E) / de of 
ionized electrons in a p -I- H2 collision are presented in 
Fig. [5^ as a function of the electron energy e and the 
impact energy of the antiprotons E. As has been men- 
tioned before, the disadvantage of the closure approach 
in Eq. (fTO)) lies in the loss of the detailed electron-energy 
information of the transitions probabilities which is of 
relevance to the electron-energy spectra (cf. |56j]). There- 
fore, the presented results may be interpreted as electron 
spectra for a fixed internuclear distance rather than in- 
cluding the full rovibrational motion of the nuclei as it is 
the case for the integrated cross sections which have been 
discussed before. The electron spectra are calculated for 
a wide electron-energy range < e < 12 a. u. and for 
different impact energies of the antiproton ranging from 
48keV to 1015 keV. The contour plot on the bottom of 
Fig. [B^ shows the corresponding level curves and gives 
therefore information on the gradient of the spectra sur- 
face. It can be seen that within the whole impact-energy 
range the electron spectra decrease smoothly and mono- 
tonically for increasing e. Considering small electron en- 
ergies e < 2 a. u., the spectra fall off strongly in view 
of the logarithmic scale for all impact energies. Within 
this e interval. Fig. [B^ shows that the smaller the im- 
pact energies E the larger the values of S{e,E). How- 
ever, for larger e this uniform trend starts to cease. For 
e > 4 a. u. the overall decrease becomes weaker. Though, 
the electron spectra for small E start to decrease again 
very strongly where the fall-off of the spectra is the steep- 
est for the smallest E. Consequently, in the intervals of e 
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Figure 6: (Color online) Electron-energy spectra surface 
S{e,E) — dcr(e, £')/de given in 10"^'' cm^ / a.u. as a func- 
tion of the electron energy e in hartree and the impact energy 
of the antiproton E in keV. a) p -I- H2; b) p -I- H2. 



and E considered here, the largest value of S{e, E) for a 
given e moves from E = 48 keV at e — to E ^ 200 keV 
at e = 12 a. u. 

Cuts ^(e) of the same electron-spectra surface for ten 
different antiproton impact energies £^ = 48, 67, 95, 132, 
186, 260, 367, 515, 723, and 1015 keV are presented in 
Fig. [71 The inset shows the same S{e) curves in an inter- 
val of small electron energies < e < 2 a. u. Thereby, the 
scaling of the y axis of the inset is kept as it is the main 
graph. The ordering of the S{e) curves in the inset is 
according to their impact energy E, i.e., the uppermost 
curve is the one for the smallest (48 keV) and the lowest 
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Figure 7: (Color online) Electron-energy spectra curves 
S{e) — da{e) / de for p + H2 as a function of the electron 
energy e at £ = 48, 67, 95, 132, 186, 260, 367, 515, 723, 
and 1015 keV. The inset shows the spectra for the range 
< e < 2 a. u. without scahng the y axis. The curves in the 
inset are order accordingly to the impact energy E. The up- 
permost curve belongs to the smallest (48 keV) and the lowest 
curve to the highest (1015 keV) impact energy E. 



curve the one for the largest (1015 keV) impact energy E. 
It can be seen that no crossings of the electron-spectra 
curves S{e) occur in this low electron-energy regime. 

In contrast to the behavior for small e shown in the 
inset the S{e) curves start to cross each other at higher 
electron energies. The curve for E — 48 keV starts to fall 
off much steeper than the other S{e) curves for e > 3 a. u. 
and therefore crosses all lower lying curves. Its first cross- 
ing takes place at e « 3.19 a. u. while its last crossing oc- 
curs at e w 6.13 a. u. with the curve for E — 1015 keV. 
The other electron-energy curves for higher antiproton 
impact energies share the same characteristics, namely, 
that the curve with the largest values of S{e) in a cer- 
tain e range starts to fall off steeper than all other lower 
lying spectra curves for higher impact energies. Though, 
with increasing impact energies E the decline of the S{e) 
curves starts at larger e and gets less steep. 

For comparison to the antiproton results in Fig. [5^ 
an electron-energy spectra surface S'(e, E) is also pre- 
sented for p + H2, i.e., for proton impact, in Fig. [5Jd. 
The electron spectra are given for the electron-energy 
range < e < 3.5 a. u. and for proton impact energies 
from 48keV to 310 keV. In general the values of S'(e, E) 
decrease for larger e. However, the most striking feature 
of Fig. [6j), in contrast to the case of antiproton impact, 
is the existence of local maxima of the spectrum curves 
S{e) for a given impact energy E which are also visible 
in the contour plot on the bottom of the figure. The 
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Figure 8: (Color online) Electron-energy spectra curve S{e) = 
da(e) / de for p -|- H2 as a function of the electron energy e at 
_E = 48keV. Theory: solid curve, present results for electron 
loss. Experiment: filled circles, results for ionization, Gealy 
et al. [5y|. The energy e of an electron with the velocity of 
the proton Vp is indicated by the vertical line. 



position of the peaks of ^(e, E) varies with the impact 
energy E. At the center of the maxima the ratios of the 
two energies E and e are such that the classical velocities 
of the proton Vp and of the electron Ve are equal, i.e.. 



/ O 7^ 



which can be reformulated as 
_ E _[ 



ivpr 



(24) 



(25) 



where Mp is the proton mass. The accuracy of this state- 
ment is demonstrated in Fig. [8] where the present 5'(e) 
spectrum curve for protons with an impact energy E = 
48keV, i.e. Vp = 1.3856 a. u., is shown as solid curve. The 
maximum of ^(e) is located at e = i (^p)^ — 0.96 a. u., 
also indicated by the vertical line. 

The occurring maxima can be explained with the sim- 
ple picture of the electron-capture process where the elec- 
tron is captured by the proton and moves basically with 
the momentum of the projectile. Therefore, the veloc- 
ity of the captured electron relative to the II2 molecule is 
given by the velocity of the projectile, namely the proton, 
as well as the electron velocity relative to the moving rest 
frame of the projectile. Since both contributions to the 
electron momentum can be oriented in different direc- 
tions the peaks of the electron spectra ^(e) are centered 
around the energy e which corresponds to a free electron 
with the velocity of the projectile, cf. Eq. ([25|) . It may 
be mentioned that the capture peaks get much less pro- 
nounced for higher impact energies. This is, first, due to 



the diminishing probabihty of capture for larger E and, 
second, due to a broader e distribution of the captured 
electrons. 

If the discussed maxima of S{e, E) in Fig. ^p are re- 
moved one is left with a smoothly decreasing electron- 
spectra surface 5nc(e, E) for increasing e which is similar 
to the one for antiproton impact in Fig. [6^. This modi- 
fied Sac{^, E) for proton impact may be interpreted as the 
electron-energy spectrum surface where the electron cap- 
ture by the projectile is excluded. In Fig. [5] the present 
S{t) curve for a proton impact energy E = 48keV is 
compared with experimental data by Gealy et al. [59| for 
which capture is excluded. The comparison shows that 
except for the e regime where capture is the dominant 
process, i.e. 0.6 > e > 1.6, the present results agree 
with the experimental data though they underestimate 
the experimental findings for high electron energies. The 
integral of the difference between the present and the ex- 
perimental curve over e ( « 2.2 x 10~^^cm^) yields the 
capture probability for £■ = 48 keV calculated by Shingal 
and Lin [131 (~ 2 x 10"i<5cm2) to a good extend. The 
reason for the structures of the theoretical curve for en- 
ergies close to the capture peak is not exactly known. 
It is likely that they originate from the finiteness of the 
numerical description. 



F. Comparison of the models Vmod and V^cai 

To the best of the authors' knowledge the only existing 
calculation for p -I- II2 collisions was performed by Ermo- 
laev who used the potential V^cai to describe the target 
which is basically an atomic hydrogen target Hgcai with 
a scaled nuclear charge Zn — 1.09. His results, shown 
in Fig. [31 are not conform with the experimental data 
and the present findings. In order to find out why the 
one model Knod yields much better results than the other 
and whether the same disagreement occurs also for p im- 
pact the same cross sections were calculated again for p 
and p collisions but using V^cai in order to describe the 
target. The resulting cross sections for p and p impact 
multiplied by the factor two are also shown in Figs. [3] to 
[S] as dash-dotted curves. In what follows three remarks 
shall be made concerning the results of the calculations 
using the scaled hydrogen potential V^cai- 

First, it is obvious that the present results for ioniza- 
tion in p + Hscai collisions shown in Fig. [4] clearly devi- 
ate from those of Ermolaev [51|. It is astonishing that 
the latter results by Ermolaev better match the present 
data for unsealed atomic hydrogen than the present and 
experimental data for molecular hydrogen. No detailed 
information is given in 51] concerning the employed ba- 
sis set and the convergence of the calculations. In very 
similar studies by Ermolaev [6l|, |62], however, a two- 
center Slater- type orbital expansion with 51 basis func- 
tions were applied to describe the collision process be- 
tween p and p projectiles and (unsealed) hydrogen atom 
targets. It may be mentioned that the quality of the con- 



tinuum description in the calculations by Ermolaev has 
been put into question by other authors [40, mM i espe- 
cially in the so-called 'polarization region', i.e. between 
~70and500keV. 

Second, the present ionization cross sections using 
Vmod and V^cai as target potentials both yield, especially 
for E > 100 keV, comparable results which can be seen 
in Figs. [31 and [H for p and p impact, respectively. De- 
viations become visible for E < 100 keV. The similar 
behavior may be explained with their comparable ioniza- 
tion potentials /mod = 0.5976 a. u. at i?„ = (i?„) and 
/p[Hscai]= 0.5945 a. u. 

Third, the present cross sections for excitation in Fig. 
\5\ differ, however, considerably for Knod and T^cai- To 
the best of the authors knowledge no literature data ex- 
ists to compare the present results with and therefore 
to judge which of both models is superior in describing 
excitation of II2 molecules. On the other hand it has 
been observed that the excitation cross sections of alkali- 
metal atoms depend considerably on the energy differ- 
ence Ae between the energetically lowest dipole-allowed 
p states and the ground states [4J|, i.e., the excitation 
energy. In this context, it shall be noted that also in the 
present investigation the dipole-allowed transitions from 
the ground state to the bound states with angular mo- 
mentum / = 1, namely the p states, play, especially for 
E > 100 keV, a dominant role. The cross sections for 
excitations into bound / — 1 states are also shown in Fig. 
[SI for Vmod as thin short-dashed curve and Vscai as thin 
dash-double-dotted curve. In [53 the excitation energies 
Ae for dipole-allowed transitions are compared for the 
two models Vscai and Vinod- Therein, it turns out that the 
excitation energies calculated with Vscai are smaller than 
those for Vmod throughout the i?„ range which is consid- 
ered here. The substantial differences of the excitation 
cross sections in Fig. [5l can therefore be understood by 
considering the diversity of the Ae curves for both model 
potentials, namely, the lower excitation energies for l^cai 
lead to larger excitation cross sections compared to those 

of v;„od- 

In order to find out how well the excitation is de- 
scribed by the employed models the excitation-energy 
curves obtained with Vmod can also be compared with 
the Ae curves of exact H2 calculations. For such a com- 
parison it has to be considered that the H2 molecule can 
be oriented arbitrarily in a collision. Therefore, transi- 
tions from the II2 ground state to states with the molec- 
ular symmetries ^E„ and ^n„ are both dipole-allowed. 
The molecules are oriented statistically 2/3 perpendicu- 
lar and 1/3 parallel to the projectile momentum. Conse- 
quently, the sum of accordingly weighted excitation en- 
ergies, namely, 2/3Ae{^Uu) + l/3Ae(il]„), should be 
compared to the excitation energy of the model Vmod for 
transitions into p states. The comparison for the whole 
Rri range considered in the present work was done in 
[55| and yielded a good agreement between the consid- 
ered excitation energies of the model and the exact II2 
molecule. Therefore, it is reasonable to assume that also 
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the present excitation cross sections calculated with the 
model Knod are superior to those calculated with Kscai- 
For high impact energies the results of the model Vinod 
may even match the excitation cross sections for exact 
H2 molecules completely. 

It shall be emphasized that Vinod depends only on one 
parameter a which is determined by the ionization po- 
tential. There are no additional parameters in order to fit 
the energies or wave functions of excited states. There- 
fore, it is remarkable that in spite of the simplicity of the 
model potential Vmod it is possible to reproduce cross sec- 
tions reasonable well for ionization and excitation of H2 
molecules in strong laser fields |5J| as well as in collisions 
with antiprotons. 



IV. CONCLUSION 

Time-dependent close-coupling calculations of ioniza- 
tion and excitation cross sections for antiproton and pro- 
ton collisions with molecular hydrogen have been per- 
fornrcd in a wide impact-energy range from 8 to 4000 keV. 
The target molecule is treated as an effective one-electron 
system using a model potential which provides the cor- 
rect ground-state ionization potential for a fixed inter- 
nuclear distance and behaves like the pure Coulomb po- 
tential of a hydrogen atom for large r. The total wave 
function is expanded in a one-center approach in eigen- 
functions of the one-electron model Hamiltonian of the 
target. The radial part of the basis functions is expanded 
in B-spline functions and the angular part in a symmetry- 
adapted sum of spherical harmonics. The collision pro- 
cess is described with the help of the classical trajectory 
approximation. 

It was found that the ionization cross sections depend 
approximately linear on i?„ in the interval 1.0 a. u. < 
Rn < 2.11a. u. The dependence of (Tion(-Rn) on i?„ di- 
minishes with higher energies. Cross sections which ac- 
count for the vibrational motion of the H2 nuclei can be 
obtained by employing closure, exploiting the linear be- 
havior of crion(i?ri); ^^^^ performing the calculations at 
Rn = {Rn) = 1.448 a. u. 

The results of the calculations for electron loss in p 
+ H2 collisions agree with experimental and theoreti- 
cal data indicating the applicability of the used method. 
The present ionization cross sections for p -I- H2 colli- 
sions agree for E > 90keV with the experiment. For 
20keV < E < 80keV the magnitude of the calculated 
(Tion is still comparable to the experimental data, though 
both curves start to have a different slope. The calcu- 
lated excitation cross sections for p -t- II2 collisions were 
found to be very similar to those for the excitation of a 
single hydrogen atom by antiproton impact. 

An electron-energy spectrum surface ^(e, E) for p -I- H2 
collisions is presented for a wide electron-energy range 
< £ < 12 a. u. and for impact energies 48 < -E < 
1015 keV. In the interval e < 3 a. u. the electron-spectrum 
curves ^(e) for fixed impact energies E are smooth curves 



which do not cross. The S{e) curves are ordered accord- 
ing to the corresponding impact energy with decreasing 
magnitude for increasing E. For higher e crossings of the 
S{e) occur. Thereby, it is always the uppermost curve 
which crosses all lower lying spectrum curves S{e) which 
belong to larger E. The present electron-energy spec- 
trum surface S{e, E) for p -I- II2 collisions also includes 
the electron-capture by the projectile which manifests it- 
self in local maxima of the spectrum curves S{e) for a 
given impact energy. The position of the peaks of S{e) is 
given by e = E/Mp = {vpf/2. 

A comparison of the used model potential V^wd with a 
scaled hydrogen atom with comparable ionization poten- 
tial yields similar ionization cross sections. Therefore, 
the ionization process appears to be mainly depending 
on the ionization potential. The cross sections for exci- 
tation, however, differ notably which may be explained 
with the differing binding energies of the dipole-allowed 
bound states in both models. Since the excitation en- 
ergies of the lowest p states of the model Vmod coincide 
with the statistically-weighted dipole-allowed excitation 
energies for the II2 molecule the model Knod is considered 
to be superior to the description with a scaled hydrogen 
atom. 

Concerning the applicability of the used model poten- 
tial V^iod it was demonstrated that it is suitable for de- 
scribing ionization in p -f II2 collisions at impact ener- 
gies E > 90keV. Furthermore, the model is capable of 
determining the dependence of the cross sections on the 
internuclear distance. Even the calculation of excitation 
cross sections seems to be meaningful. Thereby, it has to 
be emphasized that besides the one parameter a which is 
directly determined by the ionization potential no addi- 
tional parameter is included in the potential in order to 
fit the energies or wave functions to those of the correct 
electronic states. On the other hand, not all effects which 
may be of increasing importance at low impact energies 
can be described by the model. First, the influence of a 
second electron is solely incorporated as a screening, sec- 
ond, no dependence on the molecular orientation during 
the collision is allowed for and third, vibrational exci- 
tation which also includes dissociation is not considered. 
Therefore, it would be eligible to perform full calculations 
which take the molecular properties of the target as well 
as the two-electron effects, like double ionization or ion- 
ization excitation, into account. Such a theoretical effort, 
which is currently in preparation, accompanied with pre- 
cise measurements at low antiproton energies would lead 
to a better understanding of the p -I- H2 collision process 
for E < 100 keV. 
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